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Abstract — In this paper we introduce a modified 
mathematical model to represent spread of disease, which 
adheres to the structure and dynamics of social network. We 
also discuss the importance and effectiveness of social network 
based study in the prediction and control of epidemics. 

Index Terms — Mathematical epidemiology; Social network; 
Mathematical models; Core. 

I. INTRODUCTION 

First remarkable study in mathematical epidemiology is the 
study of smallpox by Daniel Bernoulli in 1760. But widely 
accepted models in mathematical epidemiology were 
proposed by public health physicians such as Sir Ross, R. A., 
W. H. Hamer, A. G. McKendric and W. O. Kermak between 
1900 and 1935. These models are called compartmental 
models because they are obtained by dividing the host 
population into different compartments such as Sunciptible 
(S), Infected (I), Recovered (R) etc. For a detailed study of 
this subject the text edited by Fred Brauer et al. [1] is a better 
choice. The compartmental models such as SIS, SIR etc. are 
oversimplified with the assumption that the host population is 
homogenous and contacts between members are same 
throughout the population and steady over time. As evident 
from the literature on mathematical study of epidemics, this 
simplification has led to much variation between theory and 
practice. 

Diseases spread among members of a society due to many 
reasons. Person to person contact is one among them, which 
is the main focus of our study. Our society is heterogeneous 
and human contact depends on the personal factors such as 
likes, dislikes, needs etc. There are some external factors such 
as disease monitoring and control mechanism, mass 
campaigning about the disease, availability of drugs and 
preventive medicines etc. In the initial stage of the spread of a 
disease people do not have proper information about the 
occurrence of the disease and in many cases authority many 
be weak in implementing right preventive measures. People 
may not be aware of the nature and consequence spread. Later 
as the disease progresses to an epidemic or endemic, more 
and more people become aware of the situation and they tend 
to reduce contact with infected people and in some cases 
society itself declares quarantine against diseased. As disease 
become stronger in infected people, they normally have to 
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take bed-rest which results in some kind social withdrawal. 
All the above mentioned factors greatly affect social contacts 
between infected and susceptible and hence control the 
spread of disease in the later stage. Many previous studies 
have shown the effect of contact network and human 
behavior in the spread of disease. One is the study about the 
spread of gonorrhoea in the United States [2], which showed 
60% all infections were caused by a small group of 
individuals that constitute only 2% of the whole population. 
In this study we try to incorporate some of these factors while 
modifying existing models. The outbreak of SARS in the year 
2003 also hints about the role of structure in mathematical 
epidemiology [3]. 

In the next section we review a simple models being used in 
epidemiology, the SIS model and discuss its characteristic 
features. We also discuss its disadvantage and weakness in 
explaining some phenomenon. Then we move to the third 
section in which some networks with different structures are 
selected and analyzed to prove the weakness of homogenous 
mixing models. A revised model is given in the fourth section 
and its effectiveness in the analysis of disease dynamics is 
discussed. 


ii. The SIS model 

For a person who studies the mathematics of epidemics, SIS 
model is a starting point. This model is obtained by dividing 
the host population ( N stands for the total size of the 
population) into two compartments; the Susceptible ( S ) and 
the Infected (7). Initially it is assumed that the whole 
population contains only susceptible people and a disease is 
newly introduced into the population. The disease spreads 
among the members when they come in contact other 
members. A third assumption is that people do not acquire 
permanent immunity against the disease. A person, who 
enters the compartment, returns to S after some time, which 
varies from disease to disease. So the SIS model is suitable to 
represent only those models which do not develop prolonged 
immunity. Let /? be the average rate at which a member 
contacts with others. So the total contact one person can make 
in the population is JSN. The probability of a random contact 
with a susceptible is S/N. Thus the number of new infections 
per individual per unit time is (J3N)(S/N) and the rate of new 
infection is obtained by multiplying this quantity by I. So the 
number of people who leave the compartment S to join I is 
fiSI. If y is the rate of recovery, yl is the number of people who 
are recovered per unit time. Thus we get the following 
differential equations which model the disease dynamics. 
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This model is built on the assumption that all members have 
same kind of contact with a specific percentage of other 
members. The above equations represent a network in which 
members have homogenous type of mixing with the rest of 
the population. Graph theoretic model of the network 
depicted above is strictly a regular graph. A regular graph is a 
graph in which all vertices have same number of degree. Next 
we proceed to discuss how variations in network structures 
affect dynamics of disease spread in contact networks. 

hi. Regular Graphs and Spread of Disease 


To simplify the analysis, we assume that initially one 


individual is infected ( 



and in turn he infects 


/V — 1 > 


others in unit time. 


If we fix one unit of time as period, during the next period, we 
have the following values for ~~ x s respectively. 


A/, =7V 
h =J3( A^Ul 
S\ =GV-U+y 



Again repeating this for the next period, we get a more 

complicated set of equations, which are functions of - J 
Generalizing the above steps, we can develop a set of 
recurrence relations, which give the values of 

S -ill* l 


^ / \) 1 ) 



on who is infected and not on the size of the population. Next 
we consider another contact network known as an 
n-dimensional hypercube graph. We can label the vertices in 

a hypercube, which contains 0 vertices using the 


tuple 


J , where & 


=0 orl 


. Two 


vertices are adjacent if and only if there is exactly one change 


in the values of 0 in the corresponding positions. In this 
network time taken to cover the entire network is same 
irrespective of the position from where we start. So this 
network is highly homogenous unlike the previous network. 
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Figure 1 . Three dimensional Hypercube. 

Disease can spread equally in all directions in this type of 
networks. The standard deviation of the degrees of vertices is 
almost zero in both. Clearly S. D. of degree of vertices is not a 
factor which decides the propagation of disease in a network. 



By substituting the values of ^ ~ v obtained in 
previous steps, we can show that all compartments are 
expressed in terms of the total population. 


Now we consider a strict binary tree structure (fig: 2) for the 
contact network. A strict binary tree is a binary tree, in which 
each vertex, except the leaves have two branches. Assume 

that a person occupying the position 0 is infected. He 
then infects a portion of individuals in his neighbourhood. 
These infected people infect the other people in their 
neighbourhood union and so on. In graph theory, an open 


7V(v) 


IS 


neighborhood of a vertex which is denoted by 
the set of all vertices which are connected to that vertex. The 


closed neighbourhood of a vertex, 7VCl ’ 3 is the union of its 
open neighbourhood and the vertex itself (i.e., 

) jn a binary tree if the infection is at 
the root, the disease has to pass through four steps to infect all 
members. On the other hand, if a leaf vertex is infected, eight 
steps are necessary to cover the entire population. Thus time 
taken by the disease to infect the whole population depends 


Another highly homogenous regular graph network is a 
complete graph network. In this network every member is 

connected to every other member. An - vertex complete 

graph contains 

degree 


(n — 1 ) 


edges and it is regular with 

. If a vertex in a network having complete 
graph structure is infected, then it can transmit disease to all 
other members in single step and can produce an epidemic 
within a short period. This character of the vertex is due to the 
fact that the closed neighbourhood of every member is the 
entire population. It motivates us to define the core of a 
network. It is the set of vertices in a network, such that the 
union of the closed neighbourhoods of the members of the set 
spans the entire network. Each member in a complete graph 
network is a core. Core is a very crucial idea in the study of 
epidemic. If an infection occurs in the core, it can definitely 
spread rapidly in one or two steps and thus it can cause a 
pandemic. 
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iv. Other Networks 

A. Star Structured Networks 

A star in graph theoretic terminology is a graph, which 
contains a vertex that is connected to all other vertices having 
degree one. So there is only one vertex having degree greater 
than one. This vertex is called hub in the terminology of 
network science. 



Figure 3. A star network 


A star structured network is a union of many stars. So it can 
be decomposed into small star sub-graphs 
having hub vertices 


' i ^ 


The core of the graph is 


hi 




. A star structured network is highly 
vulnerable to disease because it can spread over the whole 
network in one or two steps if the first infection is at all 
vertices in the core. On the other hand if only a part of the 
core is initially infected then the disease has to take many 
steps to spread over the network depending on the network 
topology. 


B. Random Networks 

Random networks are artificial networks. They do not 
possess any specific structure. So they are not representative 
of any natural networks. But they provide a platform for 
comparative study of structured networks. Random networks 
are identified with their degree distribution. They have either 
a Poisson degree distribution (if the network size is very big) 
or a Binomial distribution (if network size is very small). 


C. Scale-free Network 

Scale-free networks are the networks in which there are a 
large number of vertices of small degree and very small 
number of vertices of large degree. Thus the degree 
distribution obeys a power law. The large degree vertices are 
called hubs. So a scale-free network contains a small 
collection of hubs. The neighbourhood union of this set spans 
almost the entire network. So it behaves like a core. If disease 
affects a considerable number of members of the core, it 
quickly spreads in the networks. If it affects a set of members 
who have very law connectivity, the rate of spread may 
decrease and ultimately die out. So the question, which part 
of the population is affected by the disease is very important 
in the study. 


v. The Model 

In the light of above discussion, we propose a new model to 
represent spread of disease in a structured network. We assume 
that the process contains only infected and susceptible. Nobody 
is recovered permanently and birth and death does not occur in 
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the population. Let 7 be the infective period for a disease. 
This is the average time an infected remains infective. We 
divide the whole period of time in subintervals having 

length T and @ denote the 0 th interval. Let ® , § 

and b e faQ infected, susceptible and 

the open neighbourhood set of the infected members, 

respectively in the beginning of the 0 th interval. We assume 
that the population size is steady and /, is initial set of 

infected. Then, v () = - i Q . n(i 0 ) is the set of all 


members who are in contact with /, . Due to many factors 


such as weak contact, natural resistance developed by internal 
mechanism against disease etc., all of them may not be infected. 

So we introduce a probability factor here. Let fl be the 
probability of a contact resulting into an infection. Then the 


total number of infection is 


« i jsrc ^ ) i 


. By time 7 


all 


infected might have become susceptible. Thus after time 7 
we get, 




After time ® , 


I^=*ATS£I^\ 

.sy =/v — /vy, > i 


Generalizing the equations, we get, 





To get the total number of infections throughout the period of 
prevalence of disease, we can use the formula, 


I 

/=( ) 



and the probability of infection is roughly obtained by 




In the model given above, it assumed that the average time 

0 for transmission of disease from one group of infected to 
the set of individuals in its neighbourhood is same as the 

average rate of recovery 


If 




an infected 


recovers faster than he makes a susceptible person, infected. 

So strength of the disease decreases and if ^ <<7 l3 
ultimately the disease dies out. On the other hand, if 

T -’ disease spreads faster in the neighbourhoods and 


it grows gradually. If the inequality is very strong, then the 
disease turns out to be a pandemic. 


On the outbreak of a disease, through mass media and 
campaigning programmes, people are made aware of the 
disease. It affects the attitude of the people towards the disease 
and they resort to some self defensive mechanism. For example, 
they keep themselves away from the diseased person reducing 
their public contacts. Consequently, the associated network 
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shows changes in the structure and hence it becomes necessary 
to consider dynamic graphs instead of static one. 


vi. Further On Core 

Here we generalize the concept of core to incorporate the 
structure of a network and to make the idea more meaningful 
and applicable. The second order closed neighbourhood of a 

set v Vf<r ' i s the subset ^ and 

the third order closed neighbourhood is 

^ ^ ^ ^ ^ This concept can be 

generalized to obtain the n th order closed neighbourhood of 

the set 0 We define the 

span number of a graph with respect to a given 

subset s as the smallest number 0 such that 

If G = K " the complete graph containing n 

for all If 

and 


7sr n | s | = v . 


vertices, then 

where 


g k 


(i ,»> 


i »>) > = 1 


n > 2 , 




d (v) = n 


then 

{v} ) = 1 

and 

{vy } >12 


for 


all i. From the second example it is clear that for two subsets 

2 > I «(G, .S' | ) 


and 


and 


oc{G, S 2 ) 


may be 


different. We can prove the following results on closed 
neighbourhood and span number of subsets of 17 


Lemma - 1 : For any two subsets 

[ 1 ] 


and 


of 0 




[ 2 ] 

Proof of 1 : For any two vertices V| and 

Generalizing it for 

an arbitrary set, we get |v — 

Proof of 2: Applying 1.1 repeatedly n times, we get lemma 2. 

Theorem -2: For any two subsets 0 and H of 0 such 
that 151 ^ 50 , 


Proof: Since 


can span the whole graph in 
cx ( g , s x ) \ s ^ e p S But H contains more vertices than ® 


So, it may be possible to cover the graph in lesser number of 
steps. 


The above relations are important in the determination of 
spread of disease. Correct information about the higher order 
neighbourhoods of a set of infected is very crucial in deciding 
the policies of disease control. Vaccinating all members of the 
neighbourhood of infected is a better choice than doing it for 
the whole population, as a temporary, quick and cost effective 
measure. 


vie Importance Of Network Structure In Control Of 

Disease 

From the discussions so far, it is clear that the knowledge about 
the spread and related dynamics of disease is very important in 
all stages of control and management of disease. Spread of 
disease in the network is controlled by the network topology. 
While vaccinating against disease, vaccinating the whole 
population is costly and very difficult. But if the members of 
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the core are vaccinated first, authorities could effectively 
control the disease. Widely used models like SIS and SIR do 
not account for network structure. Heterogeneity is of least 
importance in these models and uses the assumption that each 
member is likely to contact with any other member of the 
population. Analyses of original data have shown wide 
differences of theory from reality. This has motivated us to 
develop a new model, in which contact of a person in his 
neighbourhood and structure are the main concern. We can 
modify it further incorporating the factors affecting dynamism 
of network to get much improved models. 
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